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Abstract 

We derive the full canonical formulation of the bosonic sector of 1 1-dimensional super- 
| gravity, and explicitly present the constraint algebra. We then compactify M-theory on a 

. warped product of homogeneous spaces of constant curvature, and construct a minisuper- 

space of scale factors. First classical behaviour of the minisuperspace system is analysed, and 
then a quantum theory is constructed. It turns out that there similarities with the "pre-Big 
Bang" scenario in String Theory. 

in 

! 1 Introduction 

One of the most fundamental problems in theoretical physics is the search for a quantum theory 
which would unify gravity with other interactions. Over the past 20 years, superstring theory 
emerged as a successful candidate for this role. It was later discovered that all five superstring 
theories can all be obtained as special limits of a more general eleven dimensional theory known 
as M-theory and moreover, the low energy limit of which is the eleven dimensional supergravity 
i-^ ■ [HE]- The complete formulation of M-theory is however not known yet. 

In a cosmological context, there is another approach to quantum gravity which was pioneered 
by DeWitt in [3]. Here, the gravitational action is reformulated as a constrained Hamiltonian 
system and then quantized canonically. The resulting wavefunction is sometimes referred to 
as the "wavefunction of the universe" [I], as it describes the state of the universe. Such a 
wavefunction is a function on the superspace - an infinite dimensional space of all possible 
metrics modulo the diffeomorphisms. Although this procedure of course does not give a full 
theory of quantum gravity, it does give a low energy approximation, which is enough to capture 
some quantum effects such as tunnelling [H Since the behaviour of the wavefunction in the 
full infinite-dimensional superspace is difficult to analyze, models with a reduced number of 
degrees of freedom have been considered. In these models only a finite subset of the original 
degrees of freedom are allowed to vary, while the rest are fixed, so that the wavefunction becomes 
a function on a finite-dimensional minisuperspace. In the early Universe, it is perceived that 
quantum gravity effects should become important, and so such minisuperspace models, where 
the degrees of freedom are the spatial scale factor in the Friedmann-Robertson- Walker (FRW) 
metric and possibly a scalar field, have been used to explore different quantum cosmological 
scenarios [i EJ H [8] . 

With the advent of superstring theory, the above ideas have been applied in the context of 
superstring theory (9j \10\ [TTJ H2] [13] , [H] . This time however, the starting point is the lowest 



1 



order string effective action, possibly with a dilaton potential or a cosmological constant put in. 
Hence compared to the pure gravity case, there are new degrees of freedom - the dilaton and 
any of the tensor fields that appear. The minisuperspace models studied in the quantum string 
cosmology setting now have the FRW scale factor and the dilaton field as the independent degrees 
of freedom. In particular, progress has been made in quantum string cosmological description 
of the "pre-Big Bang scenario" [10\ \TT\ [T2| 115] . In this scenario, the universe evolves from a 
weakly coupled string vacuum state to a FRW geometry through a region of large curvature. 
Classically there is the problem that the pre-Big Bang and post-Big Bang branches are separated 
by a high-curvature singularity. However, in the minisuperspace model for a spatially flat {k = 0) 
FRW space-time with a suitable dilaton potential, it is possible to find a wavefunction which 
allows tunnelling between the two classically disconnected branches and this solves the problem 
of transition between the two regimes. 

With M-theory being a good candidate to be a "theory of everything", it is interesting 
to see what the canonical quantization of the low energy effective theory can give, given the 
achievements of this approach in the pure gravity and superstring theory contexts. In Section 
2 below, we start with the bosonic action for eleven-dimensional supergravity and reformulate 
the theory as a canonical constrained Hamiltonian system. The canonical formulation of eleven- 
dimensional supergravity has been considered before in [16] and [T7], but here we explicitly give 
the constraint algebra, at least for the bosonic constraints. 

Then in Section 3, we reduce the system to a minisuperspace model. This is done by 
restricting the metric ansatz so that its spatial part is a warped product of a number of ho- 
mogeneous spaces of constant curvature, and the supergravity 4-form is also restricted so that 
only its 4-space components are allowed to be non-zero. In the case when only one of the spa- 
tial components has non-vanishing curvature and the 4-form vanishes completely, it is possible 
to solve exactly both the classical equations of motion and the corresponding equation for the 
wavefunction. 

In section 4, we consider the classical and quantum solutions in the cases of vanishing, 
negative and positive spatial curvature. Spatially flat M-theory minisuperspace models have also 
been considered in [18] and [19]. We carefully consider these systems, fixing the gauge and paying 
attention to self-adjointness of the Hamiltonian. It turns out that the positive and negative 
curvature cases can exhibit very similar behaviour to the string theory minisuperspace models 
described above with negative and positive dilaton potentials in the Hamiltonian respectively. 
However it also turns out that in the negative curvature case, the boundary condition which 
leads to a tunnelling effect, in fact also leads to lost self-adjointness of the Hamiltonian. This is 
a purely mathematical consequence, and certainly deserves further investigation to find out the 
what is the correct physical reason behind this. 

In section 5, we look at the case where the 4-form is switched on, the 3-space is flat, and one 
other spatial component is of positive curvature. 

We will be using the following conventions. The spacetime signature will be taken as 
( — h +...+) and all the curvature conventions are the same as in [20J. Greek indices fj,,u,p,... 
range from to 10, while the indices a, (3, 7, ... range from to 3. Latin indices a, b, c, ... range 
from 1 to 10. The units used are such that h = c = IGttG^ 11 ) = 1. 

2 Canonical formulation 

In this section we set up the canonical formalism for the bosonic sector of 11-dimensional super- 
gravity, with the field content being just the metric g^ u and the 3-form potential A, with field 
strength F = dA. 
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The action for the bosonic fields is 



S = J d n x(-g)% - ^Jf A*F- ^JA AF AF 

— I d ll r (—a)% I /?( n ) IWi-^P _ F Mi-Mn4 /? 7? 1 CTl 

where g = det [g^] and i^...^ = A^^^j. The 1 1-dimensional alternating tensor e^i— 
is defined by 



where j^i— -^n = — 77 is the alternating symbol. 

To decompose the metric into spatial and temporal parts, we use the following ansatz [3J: 



, -a 2 + /? a /T 13 a , 

^ " 1 ft lab ' ' ( j 



The inverse metric is given by 



where 7ac 7 bc = and a = 7 ab p b . 

Using this ansatz, we follow [20] to express canonically the gravitational action. 

Consider a hypersurface S, given by t = const. The future-pointing normal vector to this 
hypersurface is given by 



laa ) (4) 



and the corresponding covector is = (—a, 0), so hence n^n^ = — 1. 
The second fundamental form for £ is defined by 

Kfj, v = -h%hln Ka (5) 

where the semicolon denotes covariant differentiation with respect to the 1 1-dimensional metric 
g and hp is the projector onto £ defined by 
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Note that the sign in ([5|) depends on the convention used, so here we follow |20j. 
From ([5]) we have in particular 

K a b = ~n a -b- (6) 
Using the definition of (jlj), it can be shown that [20J 



Kab = 1 (/? a | b + (3 b \ a - 7 ab>0 ) (7) 

where | denotes covariant differentiation with respect to the metric j ab . Using the Gauss- 
Codazzi equation ([8]) below, the full eleven-dimensional curvature can be expressed in terms the 
intrinsic curvature of the hypersurface (that is, the curvature of the metric r y ab ) and the second 
fundamental form: 

#(ll) = #(io) _ KabK ab + K 2_ 2n ^n v R^ (8) 
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where K = j ab K ab and K ab = ^ ac ^ bd K cd . Hence the gravitational Lagrangian density C grav is 
given by 



Cgrav = (-9)^R {11) 

= (« (10) - K ah K ab + K 2 - 2n tx n v R^ 

= 075 (i?( 10 ) + K ab K ab - K 2 ^ + total derivative terms (9) 

where 7 = det i^i ab )- In the action, the full derivative terms give rise to a surface integral. We 
neglect it, since it does not affect the dynamics of the system. 
We now decompose the 3-form A^ up as 

A 0ab = B ab (10) 

A a bc = A abc (11) 

and correspondingly, 

Fated = ^d [a A bcd] (12) 

Foabc = doA abc — 3d[ a B bc j. (13) 

The F 2 term from the action ([1]) is decomposed as 

F^-^F^...^ = F abcd F abcd - AF ±bcd F ± bcd (14) 

where 

F± bcd = n^F^a = a- 1 (F obcd - (3 a F abcd ) . (15) 
Looking at the Chern-Simons term A A F A F, we have 

„Mi-Mn4 tp TP — »,oi-aio ( 19 2 R f), A 8r A (Ifi) 

+8r/ ai - ai0 (d A aia2aa ) A a4 _ asae F a7a8a9ai0 ) + total derivative term 

Again, we neglect the total derivative term, since it does not affect the equations of motion. 
Bringing together Q, (|14j) and (|16p . we thus have the total Lagrangian 

Ltot = j d 10 X (Cqrav + C 



tot — j U A, {J^grav ~r J^form) 

where 

Cgrav = r>a(R^ + K ab K ab -K 2 ) (17) 
Cform = 7 5 — -^F abcd F abcd + — -Fibcd-Fl 6cd (18) 



_ai...aio 



1 8 \ 

Y^2^o,ia2d[a 3 A a4a5af .^d^ a7 A a8a9ai0 ^ — \^Q^a\a%az) A aA a ba& F aiaia9ava j 



We see that the canonical fields in this system are a, /3 a ,7 af> which come from the gravitational 
Lagrangian, together with A abc and B ab which come from Cform. From the Lagrangian densities 
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T71) and (I18h we can now write down the canonical momenta conjugate to these variables: 

(19a) 
(19b) 
(19c) 
(19d) 



TT = 


dLtot _ 
dap 





7T a = 


OLtot 
9Pa,0 





P ab = 


OUot 

dB abfi 


= 


TT = 


OLtot 

Olabfl 


= -73 (K' 


„abc 

TT = 


dLtot 

QAabcfi 


1 1 „ 

6 ' x 



12 4'/ 7 A dld2d3 F di ... d7 . (19e) 

Expressions (fT9al) . (fl9b) and (fT9cl) are known as primary constraints [22J. This means that the 
corresponding "velocities" cannot be expressed in terms of the momenta, and are thus arbitrary. 

Now that we have the canonical momenta, we can work out the Hamiltonian for this system. 
The canonical Hamiltonian is given by 



Htot = J d 10 x (a fi ir + P afi ir a + l abfl TT ab + B abfi p ab + A abcfi TT abc - C grav - C forrtlj 
From [3], we know that the gravitational Hamiltonian H grav is given by 

H g rav = J d 10 x (a, vr + [3 afi ir a + dH + (3 aX a ) (20) 



with 



where 



H = ^ [K ah K ab -K 2 - R^) = G abcd 7r ab 7r cd - 7 ^ (10) (21a) 
x « = _2vr^ |b = -2vr^-7 ad (2 7M , c -7^)vr bc . (21b) 



Gabcd = \l » ( laclbd + ladlbc ~ gabled ) < 22 > 



is the Wheeler-DeWitt metric. 
Consider the remaining part 



Hform = y d 10 X (^B ab flp ab + ^afec^ 71 "^ ~ L~formj ■ (23) 

Due to the constraint (|19cp , nothing can be done with the first term, but in (A a bc,oTT a c — L*form) 
we have terms in A abc ^, but these can be expressed in terms of TT abc using (|19e|) . First, define 



TT abc = ^F ± abc (24) 



so that, from (jl9cH . 



TT 



abc „abc 



12 4 

Then from definition of F± abc (|15p . we have 



* abc + ^V abcdl - d7 A dld2d3 F d4 ... d7 . (25) 



Abcd,o = /3 a F abcd + 3d [b B cd] + 6«7 2 jt bcd . (26) 
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Using (|26l) to substitute A bc dfl for 7r bc d, we can write down the overall Hamiltonian in the form 



H 



tot 



d 10 x ( a, 7r + p a nvr a + B abfi p ab + aU + f) a x a + S^X' 



- ab 



where 



H 

r 

x ab 



n + —^F abcd F' 



abed 



+ 3^h abc v abc 



a -bed 



,,a _j_ 77a ~i 



(27) 

(28a) 
(28b) 



-M abc r + 



12 2 



0603... ai o 4 a a 



We see from the Hamiltonian (j27l) that the quantities a, /3 a and 5 a 6 are arbitrary, so we set the 
gauge as convenient. 

In order for the primary constraints (|19ap - (|19ep to be consistent with the equations of motion, 
the time derivatives of it, ir a and p ab must vanish. This corresponds to vanishing Poisson brackets 
of these momenta with H to t- Immediately this leads to the secondary constraints [22] 



H 

r 





0. 



(29a) 
(29b) 
(29c) 



Consequently, the Hamiltonian vanishes on the constraint surface. 

The new constraints ()29ap - (|29cp also have to consistent with the equations of motion. So 
their Poisson brackets with Htot must vanish on the constraint surface, or else there will be 
further constraints. Calculating Poisson brackets with Htot reduces to working out the pair- 
wise brackets between the quantities 7Y, x a an d X ab - The non- vanishing brackets between the 
canonical variables are: 



[TT,a'] = 5(x,x') [vr a ,/3' b ] = 5 a b S(x,x') 



P 



led 



\ed] 



ab 



5 {x, x') 



lab,** 



led 



-led) 



ab 



5 (x, a/) 



-A-abc < ^ 6 



Me}] 



abc 



5 {x, x') 



Here ' means that a quantity is evaluated at x', 5 (x, x') is the 10-dimensional delta function and 

rai...a k _ rai ra fe 
°b x .„b k ~ °b x -°b k - 

Before proceeding to the derivation of the Poisson brackets, we note that in general, the 
brackets are expressed in terms of generalized functions - <5-functions and their derivatives. So 
the technically correct way to handle them is to introduce arbitrary test functions and consider 
the action of the generalized function on them. 

The calculation can be simplified if we notice the following. For an arbitrary A ao , we have 



A 



defi 



x> ab K b d 10 x' 



A 



def,K 



labc 



Kbd l °x' 



3A 



[dej] 



So this implies that x ab is the generator of the gauge transformation 



and hence 



5\A = dA, 
5 A F = 0. 



(30) 

(31) 
(32) 
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Under this transformation, we have 



def 



4 



7T 



def 



T ,a6c 3 ...cio ( a' a' fi' A' A ' ,/ 10 ' 

12 3 ? "[ca^cscej^cr^cgcioj^ob x 



12 3 



defabcg 3 ...g 6 f) a 77 



and hence 



afec 



93—96 



0. 



(33) 



(34) 



Using (|32p and (|34|) . it immediately follows that all brackets involving x ab vanish identically, 
since relevant terms in each constraint involve only F and jr abc . 

Consider the brackets with x a now. After some index manipulation it is possible to rewrite 

X a s 

Xa = Xa + FabcdK bCd ~ A ahc ^ C ~ 3^7^. (35) 

However, \ hc is a lso a constraint and moreover all its brackets with other constraints vanish, so 
we can replace x a °y an irreducible constraint \ a given by 



Xa = Xa + F a bcdK bCd -3A abc TT bcd 4 



(36) 



It is hence enough to work out the brackets with x a - 

In pure gravity, we know from [3] that x a generates spatial translations. Hence x a ac ^ s on 
7 afe and TT ab as a Lie derivative, that is, for arbitrary £ a 



7T 



x' a ed io x' 
x' a ed w x' 



LfTT 



T mnp 



We can now work out the action of x a on A mnp and n r ' 

Abed, j x' a t' a d w x' = CF abcd + 38 [b {CA\ a]cd] ) = C^A bcd 
since A br d is a 3-form. For ir bcd we have 



(37) 



/ x' a ed w x' 



d a (7r bcd t: a ) ~ 3<9 a UK^) + 3£ [ V d l a 



r —bed 



(38) 



since TT bcd is a tensor density of weight 1. Therefore — x a generates spatial translations, and 
hence x a ac ^s a s a Lie derivative. Noting that Xb is a covector and 7i is a scalar density of 
weight 1, we immediately see that 



Xbijx a ed i \' 

n, [ x' a e a d 10 x' 



(39) 
(40) 
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Introducing new test functions a a and a, respectively, we have 



aZ' a d 10 xd 10 x' 



Hi c a )C d l0 x 



after integration by parts. This gives that these brackets vanish on the constraint surface. 

We are now only left with the bracket TL, T~C ■ From [3], we already know [H,TC'], so only 

need to work out the brackets [.F 2 ,7r 2 ] and [fr 2 , 7r 2 ] , since the other cross-terms vanish. After 
some lengthy calculations, which are given in the Appendix, we find that 



H,H' 



2x a 5, a {x,x')+x a , a S (x,x') 



(41) 



vanishes 



which is analogous to the untilded expression for \H,7i'\ in [3]. In particular, H,7i' 
on the constraint surface. 

Hence the full canonical description of the bosonic sector of 1 1-dimensional supergravity 
involves only three primary constraints (|19a|) - (|19c|) and the three corresponding secondary con- 
straints (|29ap - (l29cp . These constraints are first-class constraints - that is, their pairwise brackets 
vanish on the constraint surface, and they generate gauge transformations |22j . 

Consider now the quantization of this system. Adopting the same view as in [3], we will 
take it that any two field operators taken at the same space-time point commute. This way, the 
classical consistency conditions carry over to the quantum case without anomalies. So we can 
perform Dirac quantization [22] of the system. The constraints then become conditions on the 
wavefunction \&: 

= x a ^ = X° 6 * = (42a) 

This implies that Htot^ = 0, and hence from the Schrodinger equation, d^f/dt = 0. The 
behaviour of the wavefunction is completely determined by these constraints. 



3 Minisuperspace 

In general the wavefunction ^ is a function on the infinite-dimensional superspace which consists 
of 7 a6 (x) and A a ^ c (x) modulo diffeomorphisms and form gauge transformations. Behaviour in 
this infinite-dimensional space is difficult to describe, so it is useful to reduce the number of 
variables, by fixing some degrees of freedom. This way the infinite-dimensional superspace is 
reduced to a finite-dimensional minisuperspace. 

To reduce the number of degrees of freedom in the metric, we consider the following ansatz 
for the 1 1-dimensional spacetime metric: 

ds 2 n = -a (t) 2 dt 2 + e 2X '^dn 2 . (43) 
i=i 

Here each dtif is the metric of a maximally symmetric aj-dimensional space with radius of 
curvature ±1 or 0. Since the space-time is 1 1-dimensional, we also have a condition a\ + ...+a n = 
10. For each i, e x * is the scale factor of each spatial component. Thus the only remaining degrees 
of freedom which remain from 7 ab are the X 1 . Such an ansatz was used before in [23], among 
others, to set up a cosmological minisuperspace model, however here we take these ideas further 
to write down exact solutions of this model in certain cases, and we also consider the case with 
a non-vanishing 4-form, which is particularly relevant for M-theory. 
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For definiteness, suppose a\ = 3 and consider the following ansatz for the 4-form: 



F a /3 7 s = X°(t)e a ^ s . (44a) 
Ffivpa = otherwise (44b) 

where £ a p^& is the volume form on the 4-space with metric ds 2 = —dt 2 + dQ 2 . A similar ansatz 
has been used in [23]. With this ansatz, the degrees of freedom A abc are reduced to just X° (t). 
We use this notation to explicitly highlight the fact that this degree of freedom will also be part 
of our minisuperspace, on par with the gravitational degrees of freedom X 1 for i = 1, ...,n. 
The second fundamental form K ab is given in this case by 

Kab = -^a _1 7afe- (45) 
From the metric ansatz, we immediately get 

K ab K ab = a- 2 J2<*i( xi ) 

i=l 

K 2 = a~ 2 V 2 . 



where we have defined 

n 

V = J2"i xi - (46) 

i=l 

Hence we have 

where 7 = det ( l J a j ) ) is the determinant of the normalized spatial metric ^ ab . 

With the ansatz (|44a|) for the 4-form, the Chern-Simons term in the action (pQ) vanishes and 
the F 2 term becomes 

—F^-^F. , A = -I a - 2 e - Ml (X°Y . 
48 ^•" M4 2 V / 

If we assume spatial sections of finite volume, for simplicity we can normalize this volume to 

be unity. Thus, rewriting the action in terms of the new variables X° and X\ and integrating 

out the spatial integral, we obtain the action for the minisuperspace model S mss : 



S m ss = fdt ^ (f> (iff - V 2 + \e- 2 ^ xl (X°f^j +^e 2V R^ 



(47) 



where we have defined 

fi = ae~ v . (48) 

It can be shown explicitly that the equations of motion which are obtained from this action 
are equivalent to the equations obtained when our ansatze for the metric and the 4-form are 
substituted into the full field equations for supergravity. In particular, note that the equation 
of motion for X° is 

j t (a- 1 e v ~ 2a - xl X Q ) = 0. (49) 
The field equation for the 4-form is 



V^F^ = -±=d^{^-gF^) 



a ^e- v d (ae v F 0,/p(7 ) + -^=d a (gfF aupa ^J 
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where g\ is the determinant of the metric d£l\ for the 3-space. The second term in the sum 
vanishes due to the ansatz (|44a|) , and the remaining equation is precisely equivalent to 
From our metric ansatz (|43p . the general form of the spatial Ricci scalar i?( 10 ) is 



1=1 

where fej = ±1 or 0. With this, the action (j47j) becomes 



S mss = / rft L- 1 (j2 «< (**) 2 " ^ 2 + ^" 2aiXl (*°) 2 ) + £ ^ («* " e2(y ^ l) 

L \i=l / i=l 

(50) 

In the integrand above, we have a quadratic form in the X A for A = 0, 1, n. Let be the 
corresponding minisuperspace metric such that the Lagrangian is given by 

Lmss = li {^- 2 \g AB X a X b - (51) 

where V is the effective potential given by 

n 

i=i 

The Hamiltonian is given by 

H = n {^G AB P A P B + V\ (52) 

where the P A are conjugate momenta to the X A and G AB is the inverse metric satisfying 
G ab Gbc = 8 c- H ence the canonical form of the action is 



X A P A - pL ( l -G AB P A P B + V 



(53) 



As in the general case, the Lagrange multiplier ji enforces the Hamiltonian constraint H = 0. 
First class constraints are generators of gauge transformations, but in this case the Hamiltonian 
generates time reparametrizations, so the gauge transformations in this case are simply time 
reparametrizations. Therefore this constraint gives rise to invariance under time reparametriza- 
tion The gauge transformations generated by H are given by ^Ijg\ 

8X l = e[X\H] 
5 Pi = e[Pi,H] 
. de 
^ = 

We have a freedom of how to choose fj, and there are some natural choices for /i, but for now let 
us write \x most generally as 

H = e~ 2f (55) 

for an arbitrary function / (X A ^ [23]. The Hamiltonian equations obtained from (|53p in the 
gauge (|55|) are equivalent to the system with the Hamiltonian given by 

Hf = \ (G f ) AB PaPb + e~ 2f V (56) 
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together with the crucial constraint 

H f = 0. (57) 

where 

' G f) AB = e -2f G AB_ 



Thus effectively, different gauge choices corresponds to different conformal transformations of 
the minisuperspace metric. 

Quantization transforms the Hamiltonian constraint (I57p into the Wheeler-DeWitt equation 

®- 

H f ty f = (58) 

where is the operator corresponding to (|56p . and tyf is the wavefunction in the /-gauge. A 
general prescription for the quantized Hamiltonian operator is 



H f = F (X) (^~& f + aH f + e- 2f V^j 



where is the Laplace-Beltrami operator of the minisuperspace metric G f , K f is the Ricci 
scalar of the metric G* , a is a constant and F (X) is just some function of the X . Such an 
operator but with F = was first suggested by Hawking and Page [2?]. This operator is covariant 
under general coordinate transformations on the minisuperspace, which is the main reason for 
choosing this operator ordering when quantizing the Hamiltonian. We want the equation (j58H 
to be equivalent for any choice of /. This is true if and only if (|58p is equivalent to 

Hty = (59) 

where H = H^ =0 , ty = tyf =0 . From the theory of scalar fields in curved spacetimes [2H] we set 

$f = e bf^ (60) 

and then in order for (f58j) and (j59l) to be equivalent, it is known that we need a = and 
b = To construct a suitable Hilbert space, we need an inner product which would also be 
invariant under change of /, and also in which H* would be hermitian. It turns out that then 
the measure on the minisuperspace is given by 



dujf = e - 2f J\Gf\dX°...dX n (61) 



with \G f \ 



del [G f AB 



(62) 



. With this choice of the measure, F (X) = e 2 ^ . Hence overall, 

H f = e 2f (~& f + aK f + e~ 2f vj 
and we can say that and H are related by 

The (non-gauge fixed) inner product is then 

(ty{,ty{^ = j m{*m{duj f (63) 
and it is indeed invariant under changes of /. 
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From the inner product (|63[) we see that the momentum representation is of the form 



Pa = -i (d A + ^d A VG + ^d A f\ (64) 

where the extra term is chosen such that P A is hermitian with respect to this inner product. 

As we have seen, different gauge choices correspond to different definitions of the time pa- 
rameter in the system, so in particular there are two gauges which will be most useful for us: 

• The gauge / = 0, which corresponds to the choice [i = 1 and hence a = e . This choice 
leads to greatly simplified calculations. In [18], the time parameter in this gauge is referred 
to as gauge proper time, t g = t — to. 

• The gauge / = ^V, which corresponds to the choice /i = e~ v and hence a = 1. The time 
parameter in this gauge is the cosmic proper time, t c given by 



4 Minisuperspace solutions with a trivial 4-form 

Now consider a special case of the above scenario. Here we will consider solutions with a trivial 
4-form and we will take the spatial Ricci scalar to be 

where K 2 = 2kia\ (a\ — 1). This means that only one spatial component of the space has 
non-vanishing curvature, and the other components are flat. In particular, this special case en- 
compasses the scenario where the external 4-dimensional spacetime has a Friedmann-Robertson- 
Walker metric with k = — 1,0, +1, and the 7-dimensional internal space is a Ricci-flat compact 
manifold. This case is of interest from a cosmological point of view and also from the point of 
view of M-theory special holonomy compactifications. Hence the Lagrangian is now given as 

L f = -(:{,X'X' + ±K 2 e- 2 fe 2 ( v - xl ) 

2 y 2 

and the corresponding Hamiltonian 

H f = - {(; { f p.i'j - x -K 2 e- 2 U 2{ y- xl ). 

for i,j = 1, ...,n. In order to obtain explicit solutions in both classical and quantum cases, it is 
necessary to diagonalise the minisuperspace metric G* in such a way that V — X 1 becomes an 
independent variable. This is achieved by making the following change of variables: 

Y 1 = bi [(ai-l)X 1 + a 2 X 2 + ... + a n X n ] 

Y 2 = b 2 [( ai +a 2 -l)X 2 + a 3 X 3 + ... + a n X n ] 

Y 3 = b 3 [{ ai +a 2 + a 3 -l)X 3 + a 4 X A + ... + a n X n ] 

Y n = b n [{ai + a 2 + ... + a n -l)X n \ 
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where the coefficients bi are defined by 



-i 



Then 



and moreover 



b\ = 2ai (ai - 1) 
b\ = 2a 2 (ai - ly 1 (ai + a 2 - ly 1 

b\ = 2a n (ai + ... + a n _i - l)" 1 (ai + ... + a„ - . 
Gtx'X* = e 2/ - ( Y 1 ) " + f Y 2 ) +... + ( V ' 



b n Y r ' 



V= l - (hY 1 - b 2 Y 2 



We can now write down the Hamiltonian 

1 



H„ 



-if 



2 , 2 i , 2\ t>-2 2b7 1 Y 1 
P1+P2 + ■-+Pn) ~ K e 1 



where the pt are momenta conjugate to the Y l . The constraint H mss = becomes simply 



(- P l+p 2 + ...+pl)-K 2 e^ Y =0 
Taking into account the constraint (|68|) . the classical equations become 



e 2 ?p 1 



pi = A" 2 6i ^"Ve 2 ^" y Y 1 -- 

Pj = = e~ 2 fpj 



(66) 



(67) 



(68) 



(69) 



where j = 2, n. Since the "potential" does not depend on Y l , we get that the pi are constant. 
Also here we see that in order to be able to solve these equations easily, it is convenient to choose 
a gauge time parameter r given by 



~ 2 fdt. 



(70) 



to 



Hence in the gauge / = we have r = t g . Changing the time parameter, the equations of 
motion simplify drastically, becoming 



p 1 = K 2 bi 1 e 2b i Yl Y 1 



-Pi 



p j = 



Y 3 =Pi 



(71) 



where the dot now denotes time derivatives with respect to r. These are the same equations we 
would get with time parameter t g . Thus we get solutions for j = 2, ...,n: 

yi = PiT + y> 

where Yq are constants. Since all momenta except p\ are constant in r, we can rewrite the 
Hamiltonian constraint (|68|) as 

(72) 



K 2 e 2b i Y =e~p\ 



where £ is a constant given by 



e =pi+ ...+ P 2 n . 



Using (|72|) . the equation of motion for p\ becomes 

pi = k i {e-pi) 



(73) 
(74) 
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We have thus seen that after a change of variables on the minisuperspace, the classical 
minisuperspace system is described by equations ([72]) . ([71]) and the relation between Y 1 and 
p\. Essentially these are equations of motion of a particle moving in the potential — ^K 2 e 2bl Y 
constrained so that the total energy vanishes. Apart from the initial conditions, the solutions 
depend on the curvature parameter K 2 . In fact, from ()72p we see that the sign of K 2 affects the 
nature of equation (I74p and hence the qualitative behaviour of the solution. 

A similar system is considered in |29]-|31j. where the dynamics of scale factors is studied 
in the presence of wall potentials near a cosmological singularity, giving rise to "cosmological 
billiards" . 

From (|66|) . the volume factor e v is given by 



v 

e = exp 



l(b 1 Y 1 -b 2 Y 2 -...-b n Y n ) 



= Ae^ Y \~^ T (75) 

where A = exp (6 2 Yq + ... + b n Y£) and p s = b 2 p2 + ••• + b n p n . 

We now proceed to the quantization of the minisuperspace model. The canonical variables 
in the minisuperspace are now Y l for i = 1, ...,n, and the corresponding momenta pi for i = 
1, ...,n. The minisuperspace metric is now conformally flat, and is fully flat in the gauge / = 0. 
Eventually we will set / = 0, so we will disregard terms involving /. Hence in our general 
expression for the Hamiltonian operator (|62p . the minisuperspace curvature term 1Z vanishes, and 
the Laplace-Beltrami operator reduces to the flat wave operator, and moreover the expression 
for momentum operators ([64"]) reduces to 

Pi = -id Y i ■ 

Hence the Wheeler-DeWitt equation (|58|) for this model is simply 

(-d yl + d Y2 + ... + d Y n) f + K 2 e 2b i lYl * = 0. (76) 
This equation separates, and we get 

d 2 Y1 G kl - [K 2 e 2b ~^ Yl - kf) G kl = (77) 

where 
Here k\ is given by 

kj = k 2 + ... + k 2 n (79) 

and the ki for i > 2 are eigenvalues of the momenta pi. Note that (]77p is the precise quantum 
analogue of the classical constraint (|72|) . Moreover, it can be viewed as a one-dimensional 
Schrodinger equation with an exponential potential K 2 e 2bl Y± . 

Let us now discuss gauge fixing in this system. Consider the following change of variables: 
for j = 2, n let 

= J- Pe* = \pI (so) 

For these variables, the equations of motion become 



m = e lk2Y2 ...e lknYn G kl (Y 1 ) . (78) 
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and hence the Hamiltonian H mss in these variables is given by 



H - e~ 2f 

11 mss — c 



■ ? ?-w^ +P , + ... + ^ 



(81) 



In the reduced phase space method, we take the gauge choice 

£n - t = 0. (82) 

From the equations of motion this further imposes t = r, and hence / = 0. Hence we get the 
gauge proper time. The effective Hamiltonian is now 

Heff = ~\v\ ~ \K\^ + ^ + ... = - p(n (83) 

where the p^. > 0. Thus the gauge-fixed Hamiltonian now does not vanish in general. So when 
quantizing we precisely get the Schrodinger equation 

where 

Heff = ~ + id, 2 ... + id^ + l -K^^ (84) 

and so the solutions are hence 

$ = e ik ^...e<-^^e ik - t G kl (Y 1 ) (85) 

with G kl (Y 1 ^) satisfying ([77|) . So the gauge-fixed solutions have essentially the same form as 
([75]) . but with the condition ([82]) imposed and with ^ = ^ for % = 2, n — 1. 

Alternatively, we can use the Faddeev-Popov method. From (|8ip . the Wheeler-DeWitt 
equation is 

-\d\, + tffe... + id,^ + id, n + \K*^ = 

and the solutions are 

^ = eift.^'ti^-ie^^Gi! (Y 1 ) . (86) 
The full gauge fixed inner product is given by 

<*i|* 2 > = J dY x d£ 2 ,..dij>\ {y\z 2 ,...,Q5(@)a fp * 2 (y\ &,...,£„) 

where = is the gauge condition and App is the Faddeev-Popov determinant. For O = £ n — t, 
which gives the gauge condition (|82p . App = 1, so the gauge fixed inner product is 



(*i|*2> = J dY l d^ 2 ..AC n _ 1 ^l(Y\^,...,^ 1 ,t)^ 2 (Y l ,C 2 ,...^ n _ 1 ,t) (87) 

giving a positive definite Hilbert space. The solutions (f%5"]) and ([56]) are basically identical, and 
the gauge fixed measure derived using the Faddeev-Popov method is precisely the measure of 
the reduced space. Hence the two methods are equivalent. 

In any case, the key non-trivial part of a solution of the Wheeler-DeWitt equation is the 
function G kl {Y l ) which is a solution of equation ([77]) . which we can rewrite as an eigenvalue 
problem 

HyG kl = -A^G^. 
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where Hy is the operator 

H Y = d Yl -K 2 e 2b " Y \ 

Thus operator seems to occur in many different settings, and as such has been quite well studied. 
In particular, it appears in Liouville theory [32], [33] and also appears in models such as rolling 
tachyons [33])[35]> [36]. Interestingly, various minisuperspace models also contain the same type 
of equation, even though the potential is derived from different perspectives [10] , pT] , [T3] . Setting 
z = \K | 6ie b i Y , we get 

^ + ^-(^(K 2 )z 2 -b 2 k 2 )H = 0. (89) 

For K 2 < 0, this is Bessel's equation with an imaginary parameter, and for K 2 > this is the 
modified Bessel's equation also with an imaginary parameter. Hence the solutions of (|88p are 
linear combinations of appropriate Bessel functions. 

The operator Hy is clearly hermitian on the domain Dq of smooth functions with compact 
support, but this is not enough for a full definition of a self-adjoint operator. To construct a 
self-adjoint operator, we follow the general theory as set out in [37]. First we take the domain 
Vq of the closure of Hy on Vq, and work out the deficiency indices of Hy. In our case 
this corresponds to solving equation (|88j) for eigenvalues dbi and determining the number of 
independent square-integrable solutions in each case. General results say that there exist self- 
adjoint extensions of the operator if and only if n + = n~ , and moreover the operator is already 
self-adjoint if and only if n = 0. So to find the self-adjoint extensions we need to solve 

H Y <j) = ±i<f> (90) 

For K 2 > 0, the independent solutions of (|90f) are modified Bessel functions of first and 

3 1 

second kind - J„± (z) and K v ± (z), respectively, with r] + = e*" and rf~ = e± m . However all of 
these solutions are unbounded, and hence clearly not square-integrable. Thus in this case, the 
deficiency indices both vanish, and so the operator Hy is self-adjoint. 
For K 2 < 0, the independent solutions of ("9"""]) are 

4>f (z) = J v ± (z) and 4>f (z) = J_ v ± (z) , 

where J u (z) are Bessel functions of the first kind. In this case, only <fr^ and (f)^ are square- 
integrable, hence the deficiency indices are n + = n~ = 1, and hence there is a one-parameter 
family of self-adjoint extensions defined by 

(J) G Vq, a G cX (91) 

where 9 G (0, 1] is the parameter which defines the extension. As pointed out in [33], the above 
prescription for the self-adjoint extension domain defines the asymptotic behaviour of functions 
in T>0 since for Y 1 — > +oo, all eigenfunctions of Hy have a slower rate of decay than functions 
from Dq. Therefore the asymptotic behaviour of functions in T>q has to be compatible with the 
asymptotic behaviour of <jy[ + e 2m8 (j)2 for each 9. 

4.1 Case 1: K 2 = 

Suppose the spatial curvature fully vanishes, so that K 2 = 0. Classically this gives that p\ is 
also constant, and y 1 is given by 

Y 1 = - Pl t + Yq 1 

for a constant Yq. From ([72*]) we see that pf = £ 2 . Hence in our gauge, the solutions are rather 
trivial. 
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Note that from (I65p that the cosmic and gauge time parameters are related by 



2A 



h(bipi + ...+b n p n )t 



t c = : e 2 

bipi + ... + b n p n 

So depending on the sign of the quantity (feiPi + ... + b n p n ), t c is either always positive for all 
values of t or always negative for all values of t. Overall, this can be regarded as a generalization 
of the Kasner metric. Such solutions have been obtained many times before - both in a purely 
gravitational context [38J or as here, a special case of an M-theory model [18]. To relate to 
variables a, (5 and <f> used in [18] and |19| . set a\ = 3, 02 = 6 and 0,3 = 1, together with 



X 1 = \ a -f3-U 

2 O 



X 2 = --a- -6 
2 6 r 

X 3 = a + 2(3 + -cf). 

3 

With these relations, our Lagrange multipliers fi agree, and hence the Lagrangian (|5ip becomes 

Lmss = Z^ 1 (3a 2 - 4> 2 + 6/3 2 



1 1 



thus precisely as in ([18]) and ([H]). 

Back to our variables, in the quantum case, the gauge fixed wavefunctions which are or- 
thonormal in the gauge fixed measure are 

**„,...,*» {Y^-'tn-l't) =(2^)'^e^ yl e^2... e ^-i«n- le ^ 

where 

h\ = h\ + ... + 

4.2 Case 2: K 2 > 

Suppose K 2 > 0. From the constraint (|72p we see that we must have \p±\ < £. From (|74p . and 
using the condition on p\ , we get 



//* dp 
dt = I j = arctanh (£,~ 1 pi) 
J £, ~ Pi 



Hence 

Now Y 1 is determined by 
So 



Pl =^tanh(6^ 1 et + t ) (92) 
Y 1 = -etanh (fcr^t + to) 



y 1 = a - h log (cosh (b^Zt + to) ) 
The relation (j72[) fixes the constant ci, hence y 1 is given by 

Y 1 = -61 log (Kf 1 cosh (6^ + t ) ) (93) 

Figure Q] shows the behavior in phase space (with to = 0). We can see that this solution has 
only one branch - the negative and positive momentum sectors are smoothly connected. 
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We have the following asymptotic behaviour for Y 1 and p\\ 

t — ► +00 pi — > £ Y 1 £t = -pit . . 

t — ► -00 pi — >-£ Y 1 ~ +ft = -pit 1 j 

Let us now investigate the behaviour of the original scale factors X 1 . From definition of Y , 
X 1 = V — b± Y 1 . So let us first look at the asymptotic behaviour of V. From (|66|) . up to a 
constant we have 

V = \biY l - l -p s t. 
Thus from the asymptotic behaviour of Y\ (|94p , as t — ► ±00 we get 

V~-±t(p 8 ±hO. (95) 

Note that p 2 s — b\^ 2 < 0, but since &i£ > 0, we have p s — b±£ < and p s + b\£ > 0. Therefore, 
V — > —00 as t — > ±00, so in fact V has very similar asymptotic behaviour to Y . From (|95p . 
the behaviour of X 1 is hence easily obtained: 

X 1 ~ -i* (p, ± 6i (1 - 2bf) 



It follows that the qualitative behaviour of X 1 does actually depend on the numerical values of 
the constant momenta pj. It can easily be seen now, that all other X 1 will also be asymptotically 
proportional to t, but with different constants of proportionality which also depend on the initial 
conditions. 

By construction, the overall 11-dimensional space is Ricci-flat. However let us look at what 
happens to the intrinsic curvature from 4-dimensional point of view. The expression for the 
4-dimensional Ricci scalar is given by 



R W = -K 2 e~ 2Xl +a- 2 
2 



(01 - l)^ 1 + ai(ai - 1) (x 1 ^ 



After changing variables and the time parameter, and applying the constraint and equations of 
motion, we get 

= e - 2V Q (px) (96) 
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where Q (pi) is a quadratic expression in p\ with constant coefficients, the precise form of which 
is not important here. Since p\ is always bounded (f9"2"j) . the curvature blows up when V — ► -co, 
and as we know this does happen when t — ► ±00. So although the 11-dimensional space is flat, 
from the 4-dimensional point of view there is a curvature singularity. 

The solutions we had so far were in the gauge time parameter t. To relate it to the cosmic 
time parameter t c , we need to integrate e v . In this case 



e v 



A [K- 1 ^ sech (b^ l it + t )] ^ e^ Pst 



where as before, p s = + ••• + b n p n . For to = 0, the integral of this expression can be evaluated 
explicitly in terms of the hyper geometric function (a, b; c; z) [39J : 

t c (t) = A (2K-H) ¥l j-^—e\«*™ 2 F, (hi I - Ps ) ■ 1 + h {M - Ps) ; - e ^A 

bii-Ps v 2 4 ? 4 ? / 

From this we can at least extract asymptotic behaviour of t c as t — ► ±00 [39] , 

^ ± _ 2A(2K-y -^ e _ ¥(pa±bi0 

Ps ±b^ v ; 

where are constants, which we can choose such that Cq = 0. This behaviour is hence similar 
to the K 2 = case for p\ = ±£. We know that p s — bi£ < and p s + bi£ > 0. Therefore 
as t — ► -co, the cosmic time parameter t c approaches from above and as t — * +00, t c 
approaches c$ from below. 

Hence overall, at small t c , the overall size of the universe is very small, and the 4-dimensional 
curvature is very high, then as t c increases, the size of the universe increases and hence the 
curvature decreases up to a point, after which the universe collapses again and the curvature 
blows up within a finite time Cq". 

Now consider the quantized system. As we already know, for positive K 2 , the solutions of the 
Wheeler-DeWitt equation (|89p are modified Bessel's functions with imaginary parameter ib\k\. 
So the solutions are linear combinations of functions I^kx (z) and K^kx ( z )- If we impose the 
condition that the function be bounded, this uniquely selects (z) [39]. This choice selects 
the the wavefunction which decays as Y 1 — > +00, which is consistent with the exponential 
potential in (|77p . Boundary conditions for this type of wavefunctions have been well studied 
|12j . From the previous section, we know that the operator Hy is self-adjoint in this case, and 
as pointed out in [33J, this means that there is only one family of orthogonal eigenfunctions. 

The normalized gauge fixed stationary wavefunctions with energy E = —k 2 are 



. 1 \ 1 2ki sinh7r6iA:i ik 2 f ih i f / , h -x Y 

(2vr)~ V ^ V 



(98) 



kf — k\ + ... + 



where as before, 

From [3D], [31], we know that 

f°° —K ifM (x) K iv (x) = - [5(jj l - v )+6(t i + u)\ . (99) 

Jo x 2/xsmn-7r/i 



For definiteness we fix k n > and thus the stationary wavefunctions (|98p are orthonormal in the 
gauge fixed measure: 

{^k 2 ,...,k n \^k' 2 ,...,k' n ) = * (to - k' 2 ) ...5 (k n - k' n ) . (100) 
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These functions are not in the Hilbert space TL of square-integrable functions since they show os- 
cillatory behaviour as Y 1 — ► — oo. This problem is resolved by introducing a weight distribution 
ft (&2 ; •••) k n ), so that the function 

*i (Y\ e 2 , U-l) = j Pi K) ^ k2 ,...,k n (Y 1 ,^ Cn-l) dk 2 ...dk n 

is square integrable. From (jlOOp we see that for each i, Vl/j is indeed in H if and only if 
p i (^2, k n ) is square-integrable in the k 2 ,---,k n , and the ^ are orthonormal if p i are such. 
If we take p i to be highly peaked around (k 2 , k n ), then the qualitative behaviour of the 
corresponding ^ can be very well approximated by the non-smeared function ^?k 2 ,...,k n [42] . 

Now, a general property of K,^^ (z) is that for < z < |6iA;i|, the function is oscillatory, 
and for z > |6ifei|, the function decays with asymptotic behaviour as z — ► +oo given by 

K tblkl (z) ~ e^-^Mz-h-*. (101) 

In our case however, z = Kb\e bl Yl , so the wavefunction decays extremely fast for z > |&ifci|. 
From the classical solution (|93p . 

z = Kb ie b ^ Yl < bit, 

so the region of the minisuperspace where the wavefunction is oscillatory corresponds to the 
classically allowed region, and outside it, the wavefunction amplitude is negligibly small. 
For z — ► 0, the asymptotic behaviour of ifj&ifci (z) is [39] : 

(/ 1 \i6ifei i\ \ —ib\k\ \ 
\2 Z ) | \2 Z ) \ nn9 x 

T (ihh) + T (-ihh) ( > 



Using this, for Y 1 — ► — oo, we have the following asymptotic behaviour for ^: 

//lr/T \ibiki n jti \~ibik 



1*1 



y r(«6ifei) r(-z6ifci) y 

= + ^(+) 

where JV is a constant. Thus asymptotically, $ splits into left-moving and right-moving parts, 
\I>( - ) and $( + ) respectively, with the role of the time-like coordinate being assigned to Y 1 . These 
plane waves move along the vector (k 2 , k n -i) in the "space-like" part of the minisuperspace. 
By applying the y 1 -momentum operator p\ = —i-^rr, we find that the p\ eigenvalue for ^(~) 
is k\ and the eigenvalue for \^( + ) is —k\. The constant k\ corresponds to the classical quantity 
£ and therefore left-movers correspond to the sector of the classical solution where p\ > and 
the right movers correspond to the sector where pi < 0. Also, from (I103p . we can infer that 
fluctuates with amplitude YnhTj • 

The left-moving waves can be interpreted as reflections of the right-moving waves. Effectively, 
such a reflection is a transition from the negative momentum sector to the positive momentum 
sector. In the classical system, these sectors are smoothly connected, and in the quantum 
system, this is manifested by the fact that the reflection coefficient between the two plane waves 

is R = |* (+) | 2 / |^ (_) | 2 = 1- 

So in fact the two sectors can be regarded as reflections of each 

other at y 1 — > — oo. 

A similar behaviour was discussed in [llj, in the context of a four-dimensional gravi-dilaton 
system with a negative, specially chosen dilaton potential. Here the smooth branch transition 
arises naturally from a positive curvature background since the positive curvature term in our 
action gives rise to a negative potential in the Hamiltonian ([6j 
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4.3 Case 3: K 2 < 

Now suppose K 2 < 0. Letting K 2 = —K 2 we thus have from (I72p 

f C 2 e 2b- 1 Y 1 =p j-f (104) 

Therefore in this case we have |pi| > £, so from (|74|) . and using the condition on pi, we get 

dpi 



fef 1 fdt=f -fi^ = r 1 arccoth (rVi) 
7 J £ Pi 



Hence 

and l" 1 is given by 



Pi = £coth (b^Zt + to) 



Y 1 



-6i log 



sinh (fcf t + t ) 



(105) 



(106) 



Let to = 0. Then the phase space behaviour is shown in Figure Now we see that there are two 
branches - one for which pi > £ and t is positive, and one for which pi < — £ and t is negative. 

Y 1 





P: 



Figure 2: Phase space behaviour for K 2 < 



The asymptotic behaviour of Y , V and X 1 as t — > ±oo is the same as in the K 2 > 
case, and similarly, blows up when t — > ±oo. To give an explicit relation between time 
parameters t and t c , we first need 



e v = A 



# -1 £csch (6^t + t ) 



2°1 



e 2 



In this case, at least for to = 0, the integral can be explicitly evaluated in terms of the hyperge- 
ometric function 2 Pi (a,b;c,z). Hence have |39] 



tc (t) 



Re 2A(-2K- 1 £, 



for t > 



&i£-p s 



itf MMp) . j + 6i(6if- Pa ) . e 2 6r i C A for i < o 



Using asymptotic behavior of 2 Pi (a, 6; 6 + 1, 0) for 2 — ► 00, as t — > ±00 we get same behavior 
as in the positive curvature case [39J: 



t c 



2A [2K- 1 £ 



2°1 



Ps ± bii 



-e 2 



§t(6ip s ±£) +c ±^ 
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where are constants, and again we can choose Cg = 0. As in the positive curvature case, 
Ps — b\£ < and p s + > 0. Therefore as t — ► — oo, t c approaches from above and 
as t — ► +oo, t c approaches Cg" from below. However we also have this time that as t — > 0, 
|£ c | — ► oo. This implies that Cg must actually be negative, and the overall behaviour is that as 
t goes from — oo to 0, t c goes from to +oo, and as t goes from to +oo, t c goes from — oo to 
t£. Thus, unlike the K 2 > case, t c is unbounded. 

Hence in this scenario, for negative t c the universe collapses as t c — > Cq and for positive t c 
it expands. Moreover at t c = and t c = Cg , the 4-dimensional curvature blows up. Note that 
the scale factor X n is proportional to Y n and is hence proportional to t. But as t c — ► ±oo, 
t — ► 0, thus the volume of this component of the internal space is stabilised as t c — > ±oo. 

Consider the quantum system now. For K 2 < 0, the solutions of equation (|89|) are lin- 
ear combinations of Bessel functions - J^ib^ ( z ) with an imaginary parameter rti&i&i, where 
z = Kb\e bl yl . We know that in this case there is a family of self-adjoint extensions of Hy 
parametrized by 9 G (0,1]. Correspondingly there is a set of orthonormal eigenfunctions Xk 
parametrized by some 6 £ (0, 1] such that there is a one-to-one mapping from 8 to 6 [33J.[33j. 
These eigenfunctions are given by 

(~a\ I Ft / sinn 5 71 " (bih ~ ) \ 

x % {z) = X i^hlhk ^ + 4 J *** {z) ■ (107) 

V 2smhvr6ifei ^ sinh Itt (feifci + 2i0J j 

Note that for the special case 9 = I and 9 = 1, up to normalization, we obtain the functions 
J-ib- L k- L i Jib-^kxi which when expressed in terms of Hankel functions, are the orthonormal sets 
used in [25] and [H], in particular. 

Each of these orthonormal sets is associated with a particular self-adjoint extension of the 
Hamiltonian, so once we fix the domain of the operator, we can only use one particular set 
of orthonormal eigenfunctions. Again, these functions are not square-integrable, and as such 
strictly speaking, do not belong to H., so as before, to make precise sense of them we need to 
smear them with a peaked weight distribution p (&2, k n ), and only then we can say that they 
belong to a self-adjoint extension of the Hamiltonian. 

So overall, the normalized stationary wavefunctions with energy eigenvalue E = —k 2 are 



*2Lfc. (^2,-^-1) = -^ r e^...e«£-i€*-ix2 ) (W^ 1 ) . (108) 



As we briefly mentioned before, solutions can also be written in terms of Hankel functions 

nfcx (*) an 
as following: 



Hk ki (z) and Hf b . (z). Hankel functions are combinations of Bessel functions and are defined 



Hl(z) = ^f*) J -*(*) (109a) 
smh (utt) 

H 2 (z) = J ~ iu ^ - e ~ PnJlu Q (i09b) 

zAZ) sinh(^vr) ' 1 UyDj 



Consider the limit as z — > +oc (so that Y 1 — > +00). Then 



-1 / Z ^1 fc 71 " ■ TV 

H} blk {z) ~ J-e-e-'ie" (110) 



2 — 6l fc?T 
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If we impose the the so-called tunneling boundary condition [TT], we select only left-moving 
waves at large z so that such a solution can be written as 



VP [k 2 , .., k n ] = e^ Y2 ...e ik - Yn H} biki (Khe b " Yl ) . (112) 



Here the behaviour of the wavefunction is such that for negative Y , is mostly oscillatory, 
while for positive Y 1 , the wavefunction decays as e~ Y , which is much slower than the decay 
in the K 2 > case. In the current case, all values of Y 1 are allowed classically, whereas in the 
K 2 > case, Y 1 is bounded. This explains the different decay rates. 
For z — ► 0, the asymptotic behaviour of Hf u is [39 1: 



"i^iH Kr(f+i3 " rfit^) J (113) 



and therefore, for v = ib\k\, we get 

/ i\ \ib\k\ n \-ibiki* 

Til I \ 1 [ felfclTT \ 2 Z ) I \2 Z ) 

iblkl{ ' sinh(6iA;i7r) I r (i&ifci) T (-i&ifci) 



So the asymptotic behaviour as Y 1 — ► — oo is given by 

/ 



^ = Ne ik2Y2 ...e ik " Y " 



$(-) + *(+) 



. ibi&i / - \ — iftifei 

e b lfc ^ v 2 ; e i fcl yi + v 2 , e _ tklY . i (n4) 



V 



T (ibik\) T (—ib\ki 



where N is a normalization constant. Interpreting Y 1 as the timelike coordinate, the wave- 
function is decomposed into left and right moving waves along the vector (k,2,---,k n ) in the 
"space-like" part of the superspace. Note that k\ is proportional to the magnitude of this 
vector. 

So 



N 



r(i6i&i) 



gfeifciTT amplitude 



A? 



r(ifeifci) 



From (|114p . we get that \^\ oscillates around 

although the amplitude of oscillations is the same as for the case K 2 > 0, the fluctuation relative 
to the value of is very small for large b\k\. In this case the \?( _ ) term dominates, and \^\ is 
almost constant as Y 1 — > — oo. 

As in the K 2 > case the left moving waves correspond to the classical positive momentum, 
positive t branch, and can be interpreted as being incident from the right. The right moving 
waves correspond to the classical negative momentum, negative t branch, and can be interpreted 
as a reflection of the incident wave. The ratio of the reflected and incident amplitudes is 

M+)| 2 

C>, — J !_ — p -2bifeivr 

kl ~ |*(-)| 2 

and this gives the transition probability from the positive momentum branch to the negative 
momentum branch. But as we have seen, positive t corresponds to negative t c and vice versa. 
So we have a transition from classically disconnected negative time branch to the positive time 
branch. Thus there is finite probability of a transition between the two branches which exhibit 
very different behaviour. This is similar to results obtained in [TT] in a string theory context 
with a positive dilaton potential in the Hamiltonian. Here were obtain similar behaviour, but 
the potential naturally comes from the spatial curvature. By choosing the boundary conditions 
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as we did, we made sure that the transition is in the correct direction when compared with 
classical solutions. 

Although formally we can write down such a solution which exhibits tunnelling behaviour, 
it does not mean that it fully makes sense mathematically. Indeed, if we look at the expression 
of Hankel functions in terms of Bessel functions (1109|) , we see that in order to construct such a 
function from orthonormal functions xi , we would need to use Xk ^ or a ^ l eas ^ two different 
values of v. However these functions would lie in different self-adjoint extensions of the Hamil- 
tonian, and thus the resulting solution (|112p does not lie in any domain where the Hamiltonian 
is self-adjoint. Therefore, on the space where (|112j) belongs, the Hamiltonian is not self-adjoint, 
and hence in particular energy is not observable. Moreover, from Stone's theorem [43], we know 
that quantum dynamics if unitary if and only if the infinitesimal generator - the Hamiltonian - 
is self-adjoint. Hence in this case, with a non-self adjoint Hamiltonian we also lose unitarity of 
the system. This is certainly something in need of further investigation, because it is currently 
not clear what is the precise physical explanation for this. 

5 Minisuperspace solutions for non-trivial 4-form 

Now let us consider the case with the non-trivial 4-form. In order to be able to get solutions 
explicitly, we fix the spatial Ricci scalar to be 



where K 2 = 2k n a n (a n — 1). This means that the spatial 3-manifold which is part of the external 
4-dimensional spacetime is necessarily flat. The Lagrangian is now given as 



for A, B = 0, 1, n. As before, we need to diagonalise the metric, but unlike the previous case, 
here we need V — X n and X 1 to be independent variables. Notice that in the previous section 
Y n is proportional to X n . So if in the definitions for y 1 and Y n we replace X 1 with X n and vice 
versa, and a\ with a n , and vice versa, we get variables which perfectly fit our needs. The other 
variables can obviously remain unchanged, but we relabel them for convenience. Therefore, 
overall we get the following set of variables: 



Z° = X° 

Z 1 = ci (ai + a 2 + ... + a n - 

Z 2 = c 2 [cnX 1 + (a 2 + a 3 + ... + a n -l)X 2 ] 

Z n ~ 2 = c n _ 2 [a 1 X 1 + a 2 X 2 + ... + a n _ 3 X n - 3 + (a n _2 + a n _i + a n - I) X n ~ 2 ] 

Z n - 1 = c n _i [(a n + a n _ 1 -l)X n - 1 + a n - 2 X n - 2 + ... + a 1 X 1 ] 

Z n = c n [(a n -l)X n + a n ^ 1 X n - 1 + ... + a 1 X 1 ]. 



fl(io) = l _ K 2 
2 



-2X 



L f = L -G f AB X A X B + t K 2 e -2f e 2(v-xn 



and the corresponding Hamiltonian 
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where the coefficients cf are defined by 

c 2 = 2ai (a 2 + a 3 + ... + a n - (ai + a 2 + ... + a n - 
c 2 = 2a 2 (03 + 04 + ••• + a n - l)" 1 (a 2 + a 3 + ... + a n - 

c^_! = 2a„_i (a n - (a n _i + a n - 
c n = 2a n (a n -l)~ 1 . 

With these variables we get 

G f AB X A X B = e 2 f (e-\^ zl (z Q )\ (z 1 )' + ... + (z^ - (z")') (115) 
and moreover 

V = -\ {ciZ 1 + ... + c n . x Z n ~ x - c n Z n ) (116) 
Hence we can write down the Hamiltonian: 

U „-2/ f^ic^Z 1 ^ ,2,2, ,2 2 T^2 2c^Z n \ 

timss = ~e J [e* 1 7T + 7T! + 7T 2 + ... + 7T n _ 1 - 7T n - A e n J (117) 

where tta are the momenta conjugate to Z A . With this the Hamiltonian constraint is 

Ac~ 1 Z 1 2,2 ,2 , , 2 2 v 22c~ x Z n n / n8 \ 

e3 1 7Tq + TT 1 + 7T 2 + ... + 7T n _ 1 — 7T„ — K e = (J. (H°J 

Taking into account the above constraint (|118p . the classical equations of motion are 

tt = Z° = eW zl 



(119) 



TTj = Z 3 = 7Tj 



where j = 2, ra— 1 and the dots denote derivatives with respect to parameter r, given by ([7Up . 
For j = 2, n we immediately write down 

z j =^r+zi 

where the Zq are constants. Similarly as before, we can rewrite the Hamiltonian constraint (|118[) 

as 

eK^vr 2 + vr? + C 2 - vr 2 - K 2 e 2c " zn = 0. (120) 
where £ is a constant given by 



C 2 = tt| + - + (121) 



From the equation of motion, we have 



d-Ki 2 _ x 



d-K n 

Integrating, we get 



2ir w - -i^ cT ' z ' 

OZi n 



TrfcK 1 * = A 2 -vr 2 (122) 

K 2 e 2c~ 1 Z n = X 2 n -TT 2 n (123) 



25 



substituting these expressions into the constraint, we get 

A 2 = A 2 - A 2 . (124) 
Moreover, from equations of motion (jll9[) we have 

7T! = -^{Xj-irj) (125) 

* n = c-^A 2 -^ 2 ). (126) 

Note that from (11220 . we must have A 2 > and from (I123p also has to be non-negative. Moreover, 
from (|122p . we can induce that | vri | < Ai and similarly, |7r n | < A ra for K 2 > 0, and |7r n | > A n for 
K 2 < 0. 

Consider what happens to the curvature of the 11-dimensional spacetime. From Einstein's 
equation, we get 

= -\e- 2V v 2 (127) 

hence the curvature blows up as the volume of the space tends to zero. 

We now proceed to the quantization of the minisuperspace model. The canonical variables 
in the minisuperspace are now Z A for A = 0,l,...,n, and the corresponding momenta tta 
for A = 0,1,..., n. Unlike the case with a vanishing 4-form, even in the gauge / = 0, the 
minisuperspace metric Gab is not flat, as it is given by 

Gab = diag (e _ K lz \ 1, 1, -1 

so in particular in our prescription for the Hamiltonian operator (j62H both the Laplace-Beltrami 
operator A and the minisuperspace Ricci scalar 1Z are non-trivial. We indeed have 

1 d ( rr^^AB d 



and since the only non-zero components of IZab are 

ftoo = gC- 2 Tin = g c i e 3 1 



% = G ab Uab = Icf. (129) 



we get 

So the Wheeler-DeWitt equation is given by 

(eK^Sfo + d 2 zl + ... + - d\ n - \<?dz\ * + [K 2 e 2 ^ zn - ^aq 2 ) * = 0. (130) 

Due to the presence the 4-form term, and the resulting different choice of operator ordering in 
the Hamiltonian, this equation differs from equation (|76p in the previous section by new extra 
terms. The covariant measure on the minisuperspace is given by 

doj = dZ°...dZ n e-^ lzl (131) 
and hence the momentum operators are now given by 

7T! = -i(d z i- l - Cl ^j (132) 
tta = -id A for A + 1. (133) 
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The equation (I130p can be solved by separation of variables, so we use the following ansatz 

for 

$ = e iK z e iK 2 z* _ e i*n-iZ n - J - F [z x ) G [z n ) . (134) 
With this, the Wheeler-DeWitt equation (|130p becomes 



where 



K 2 — k 2 + ••• + ft^—i + g ac i 2 - 



Separating the variables, we get the following equations for F and G: 

2 * -la . \ Z7 / r ...2„icr 1 z 1 ,2 



^i-gCr^lF-^e^i ^ -KfjF = (135) 

d 2 Zn G-(K 2 e 2c " zn -$)G = (136) 
where K\ and n n are constants such that 

K \- K \ = k 2 . (137) 



In order to transform (|135p to the same form as (|136p . set F = e aZ F for a constant cr. Turns 

d 2 zl F-(K 2 e^ lzl -k 2 )P = (138) 



out that for a = 1 , we obtain 



where k 2 = k 2 — ^-c^ 2 . Thus unsurprisingly, the operator iJy appears again. In order to 
get bounded solutions of (I138p . we need k\ > 0, and the solutions are then (z n ) for 

z n = Kc n e c " . Now k 2 > and hence we also have k 2 > 0. Then, as we know, solutions 
of (|136p are Bessel functions J± iCnKn (z\) or modified Bessel functions Ki CnKn (zi), for K 2 < 
and K 2 > respectively, where z\ = 3ciKoe3 c i zl . As previously, for K 2 < 0, we get a 
one-parameter family of self-adjoint extensions (|9ip of the Hamiltonian. 

Similarly as before, let us discuss gauge fixing. We change variables so that for j = 2, n — 1 
we have 

1 

Ci = - % = ~ 2 A (139) 
and hence the Hamiltonian H mss in these variables is given by 



H ~ -e- 2f 

1±, mss — 2 



es i 7r + 7r 1 - 7r n - A e n + ^2 + ••• + ^Cn-l- V 140 ) 



In the reduced phase space method, we take the gauge choice £ n _i — t = 0. From the equations 
of motion this further imposes t = r, and hence / = 0. Hence we get the gauge proper time. 
The effective Hamiltonian is now 



H, 



eff 



eK^Vg + ^ _ ^2 _ ^2 e 2c„ l Z»l + n ^ + + ?r ^_ 2 = _ (U1) 



where the ir^ > 0. When quantizing, we take the part of the effective Hamiltonian which 
depends on Z°, Z 1 and Z n to be the same as in (I130p . so that overall, we get 

Heff = ~\ (eK 1 * 1 ^ + d 2 zl - d 2 n ) + \c^d z , - tffc... - id Cn _ 2 - \ (K\ 2 ^ - \ac-^ 

(142) 
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and so the solutions of the corresponding Schrddinger equation with E = hence 

with F (y 1 ) satisfying (fl36l) and G (Z n ) satisfying (fl36l) . A gain, the form of the solution is same 
as (11341) but with the gauge condition ^ n _ 1 — t = imposed and with Q = ^- for i = 0, 2, n — 1. 

Alternatively, we can use the Faddeev-Popov method. From ([8"T]) . the Wheeler-DeWitt 
equation is 



-aci 



is 



and the solutions are 

^ = e ^° e ^lC 2 ... e ^-iC„-i e | C r 1 2 1 ^( Z i) G(Z n ). (144) 
The full gauge fixed inner product is given by 

(*i|* 2 > = J du&l (Z°, Z\ Z n , ( 2 , C„_i) 5 (G) A FP * 2 (Z°, Z 1 , Z n , ( 2 , C„_ a ) 
where = is the gauge condition, App is the Faddeev-Popov determinant and the measure is 

(£j c = dZ°dZ 1 dZ n d( 2 ...d( n _ 2 d( n _ 1 e-^ lzl . 
For = ( n — t, which gives the gauge condition (|82l) . App = 1, so the gauge fixed inner product 

<*l|*2> = J du^l(Z ,Z 1 ,Z n ,( 2 ,.-,( n - 2 ,t)y 2 {Z°,Z 1 ,Z n ,( 2 ,.-,(n-2,t) (145) 
giving a positive definite Hilbert space with the measure 

dw c = dZ°dZ 1 dZ n d( 2 ...d( n „ 2 e-^ lzl - 

Again we see that the two methods give equivalent results. 

The equations we got here are all very similar to the equations encountered in the previous 
section, so we can write down the solutions straight away. The classical equation for 7ri (|125p 
gives 

ttj = -Ax tanh QcT 4 Ai* + to) (146) 
and from the equation of motion for Z 1 (??) and the constraint (|122l) . we get the solution for 

z r . 

Z 1 = -3d log (x^tto cosh Qc u \t + 1 ) V (147) 

This is very similar to the solutions for Y\ considered in the previous section for K 2 > 0. In 
particular, in the phase space this solution has a single branch. 

The solutions for Z n are exactly the same as the solutions for Y\ in the previous section. 
Thus for K 2 > 0, we have 

vr n = A n tanh (c^Xnt + *i) (148) 

and 

Z n = -c n log (KX' 1 cosh [c- l X n t + t x )). (149) 
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Similarly, for K 2 < 0, we get 

Tr n = C n coth(c~ 1 X n t + ti) (150) 
Z n = -cn log (\k X' 1 sink (c^Xnt + h)^ (151) 

where K 2 = —K 2 . In both cases, the asymptotic behaviour as t — ► ±oo is 

Z 1 ~ TXit 
Z n ~ =FA„t. 

We know that the volume parameter V is given by 

V = -[cnZ n -c x Z l -TT a t) + const (152) 



where 

so as t — ► ±oo, 



7T S = C 2 7r 2 + ... + Cn-\K n -\. 



V ~ ±^t (ciAi - c n X n =F vr s ) . 



Noting that X 1 = ^Z 1 and X n = V — c n 1 Z n } we get the asymptotic behaviour of the original 
variables X 1 and X n : 

X 1 - T-Ait 
3 

X n ~ ±i* (c ul - ^ (1 - 2c- 1 ) A n - 7T S ) . 

Thus as t — > ±oo, X 1 — ► — oo, and the qualitative behaviour of V depends on the sign of 
c± = ciXi — CnXn^fTTg. Both the 4-dimensional and 1 1-dimensional curvatures are asymptotically 
proportional to e~ 2V , so the sign of c± affects the behaviour of the curvature. Consider the 
following example. If n = 2, then the internal space is 7-dimensional, and moreover tt s = and 
A 2 = Ai. This immediately gives us that c± < 0. Hence as t — ► ±00, V — > —00 and from 
(|127p this implies that R n blows up whenever t — > ±00. 

Now look at the solutions of the Wheeler-DeWitt equation for this system (j!30p . Again, 
let us look at cases of K 2 > and K 2 < separately. In the case K 2 > 0, the normalized 
gauge-fixed stationary wavefunctions with energy eigenvalue E = —K? n _^ are given by 

= N^^^.^-^eW^K^ (Sc lKo e^ lzl ) K iCnKn [Kc n e c " zn 

(153) 

where 

2K\K n sinh 2>ttc\Ki sinh7r&i£; n 



N, 



(2tt)~tt 2 

Similarly as discussed in the case of vanishing 4-form, these functions form an orthonormal basis 
in the inner product (|145[) . 

For K 2 < 0, the normalized wavefunctions are given by 



(154) 

ic 

get a set of normalized wavefunctions for each self-adjoint extension of the Hamiltonian. 



(9) 

where XkJ are the orthonormal functions given by (|lU7p . Exactly as in the previous section, we 
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For K 2 < 0, if we take the Z n solution to be a Hankel function, then as in the trivial 4-form 
case, for Z n — ► — oo, we could decompose the wavefunction into plane waves similarly as in the 
case of the trivial 4-form. Then we would get a non-trivial reflection probability R^ n — ^ CriKjl7r 
from the right-moving wave to the left-moving wave, which would correspond to a transition 
from the 7r n > branch to the ir n < branch. However, as in previous situations, such solutions 
would still belong to a domain where the Hamiltonian is not self-adjoint. 



6 Concluding remarks 

We have first derived the canonical formulation of the bosonic sector of eleven dimensional 
supergravity, together with the complete constraint algebra. The brackets of the secondary 
constraints vanish on the constraint surface, so all constraints are first-class and there are no new 
tertiary constraints. When passing to the quantum system, the constraints become conditions 
on the wavefunction which govern its behaviour. 

By introducing particular ansatze for the metric and the 4-form we reduced the system to a 
minisuperspace model with a finite number of degrees of freedom. In a special case where only 
one spatial component has non-vanishing curvature, both the classical and quantum equations 
can be solved exactly. In the positive curvature case, whether with or without the 4-form, there 
is only one branch of the classical solution, where the universe first expands after starting out 
from zero size, reaches a maximum size, and then collapses again within a finite time. When 
the universe becomes small, the wavefunction can be written in terms of plane waves travelling 
in opposite directions. These waves can be interpreted as being reflections of one another, but 
since their coefficients are equal, the transition probability is 1. A similar scenario is considered 
in [11], but the effect that there is only one classical branch of the solution is achieved there by 
having a negative dilaton potential in the Hamiltonian, which is hard to motivate in a realistic 
superstring theory context. 

In the negative curvature case, the classical solutions give two disconnected branches, one 
of which is collapsing universe, and the other branch is an expanding universe. From the 4- 
dimensional point of view, there is a curvature singularity between the two branches. It is 
possible to choose boundary conditions such that at +oo the solution can be written as a single 
wave, but at — oo, it splits into two plane waves going in opposite direction. This yields a non- 
trivial transition probability between the branches. However, such a solution does not lie in a 
domain where the Hamiltonian is self-adjoint, and hence by choosing such boundary conditions, 
we lose the self-adjointness of the Hamiltonian. It would be interesting to investigate further the 
physical reasons for this lost self-adjointness, especially since the operator which appears here is 
mathematically the same as in the String Theory minisuperspace models, so the same problem 
should arise in those settings as well. 

Apart from the self-adjointness problems, we have seen that the curvature term and the 4- 
form term in our minisuperspace models, in terms of determining the behaviour of the solution, 
plays the same role as the dilaton potential in the gravi-dilaton systems derived from string 
theory. This is quite remarkable because these terms naturally from the supergravity action, 
whereas the dilaton potentials are put in by hand. It would be interesting to investigate what 
happens when there is more than one spatial curvature term. Such an ansatz would be a 
generalization of the Freund- Rubin solution of M-theory [2] ; where the eleven-dimensional space 
is of the form AdS^ x S 7 with particular scale factors for each component. In particular, if the 
3-space is curved, then there could possibly be more interaction 4-form term and the curvature 
term. Also, in further work, a less restrictive metric ansatz with a non-trivial moduli space could 
be studied, to see how the moduli space parameters evolve and what is the behaviour of their 
wavefunctions. In particular, it would be interesting to study compactifications on manifolds of 
special holonomy with time-dependent moduli. This could either involve compactifications on 
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general GVholonomy manifolds or maybe on a Calabi-Yau space times a circle. In the latter 
case, it could be investigated how mirror symmetry is manifested from the point of view of 
a minisuperspace quantization. 

Study of M-theory minisuperspace models seems to be a promising area where there is still 
much left to be uncovered, and which will hopefully aid us in the quest to further understand 
M-theory. 
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A Appendix 



In this Appendix we give the details of the calculations involved in deriving the expression (jJT 



for the Poisson bracket 
Expanding, we have 



To do this, we first need to know the bracket [TT abc , Tr' e ^ 9 ] 



~abc -lefq 
7T , 7T J a 



abc i „abcdi...d 7 j\ , , , f) r , 4, , , , „/e/q , u * efgai...a7 A> f)' a' 



QO 

abc * efg ai ...ar A f Qi A 

' 12 4 ' 



A dld2d3 d[ di A dBded7 ^ir e}9 + Y^?? 
32 



0i02i3^[a4 asaaa 7 \ 



+ 



12 4 



abcd 1 ...d 7 a 1 A, , , , Tr' e f3 



32 

32 

H tV 

12 4 ' 



e/gai...a 7 ( \ l r 
I \^ ai a 2 a 3 



a 5 a 6 a 7 d,a' 4 (x,x) + 5 ai a 2 a 3 ^ ( x i x ) ^[a A A a b a s a 7 



abcd\...d 7 



32 



Md 2 d 3 S ef9 d5d6dT S,d4 + §et9 d ld2 d 3 6 i x > x ') d [d4 A d 5 d 6 d 7 



w rf**f**"~ [A' aia2a3 5 K (x,x')+d(x,x')d' [a A' 

+ A aia2a3 5 ja4 (x, X ) + 5 (x, X ) d[ ai A a2a3a4 ]) 



a 2 a 3 a A \ 



32 



_ v abcef gai ...a 4 [a'^S ^ (x, x>) + A aia2a3 S, a4 (x,x>)+26 (x,x') d [ai A 



a 2 a 3 <H\ 



(155) 



So, for the 7i,H' bracket, we have: 



n,n' 



[H,H'] + 7 l 7 '-lF a6cd 7r / e/5 \A bcd ,^} -^ 7 ^F' abcd K efg \A' bcd ,n^ 



+367 27' 2^fafe C 7r; 



~abc ~te fq 



[H, H'] + 7^7'"^ ae/ ^' e/ /,a (x, x') - i^T^F'^H^d^ (x, x') (156) 



e/g 



X, X 



In the first line cross terms involving TC vanish because the form terms involve no derivatives 
of 7 a& and 7i does not involve any derivatives of 7r ab . Note that the undifferentiated <5-function 
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term in TTabc^' e f g [TT abc ,TT ,e ^ 9 ] vanishes, because n abce ^ 9 7r a b c TTefg = 0. Let ^ and £ 2 be arbitrary 
test functions. Then 



H, W'J Zx&dxdx' = J J [H,H'] Zt&dxdx' (157) 

77 , ~ 1 )^ ai --- a4 *a 2 a 3 a 4 <U O ~ W^Y F'"^ a2a3a J^ (x,x')) ^ 2 dxdx' 
+ (77 , )" l ^o# , e/ ^° 6Ce/SOl - a4 (<a 2 a 3 ^ +A™<U (x,x'))^' 2 ^' 

We know from [3] that 

y y [H,H'] i^dxdx' = j X a - ^ (158) 

The second line in (|157p becomes 

( 7 § 7 '-§i^/%' e/ff £ a (x,x') - ih^F' ae fH efg 5 <a/ {x,x')) S&dxdx? 

~ J J (^Hfg (t^^i) a & - -y-^efg (V^'^) ^ ^) 5 (x, X>) dxdx' 



Now look at the third line in (|157|) . After integrating by parts, and integrating out the <5-function 

we get 

y ^V^^^-^abciTefgAa^ (^ V4 - Z lm Z 2 ) dx = 

again because rj abce ^ 9 TT a b c TT e f g = 0. 
Thus 



H,U' 



W 2 dxdx' = J (x a + F a efg v ef9 ) (^e 2) a " dx 
= j T 2,a ~ f l,af 2) dx 



(159) 



Correspondingly, 



H,H' 



2x a S,a {x,x')+ X a , a S (x,x') 



which is completely analogous to the untilded expression. In particular, 
constraint surface. 



n,H' 



vanishes on the 
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